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Abstract 

We consider the thin layer quantization with use of only the most 
elementary notions of differential geometry. We consider this method 
in higher dimensions and get an explicit formula for quantum poten- 
tial. For codimension 1 surfaces the quantum potential is presented in 
terms of principal curvatures, and equivalence with Prokhorov quan- 
tization method is proved. It is shown that, in contrast with original 
da Costa method, Prokhorov quantization can be generalized directly 
to higher codimensions. 
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1. Introduction 



We consider free particle motion on a curved surface. Probably, 
the first quantum theory of it was proposed by Podolsky in [1]; he 
postulated that the Hamiltonian operator is H = —^Alb with A^s 
being the Laplace-Beltrami operator on the surface. One may try to 
get the Podolsky theory by some quantization procedure. Different 
methods of quantization yield various results which usually have the 
following general form: 

H = -^ALB + V,ix), (1) 

the function Vq{x) is commonly called "quantum potential". The 
Dirac canonical quantization [2] and the abelian conversion method 
[3,4] were discussed in our previous work [5] for surfaces of codimen- 
sion 1 ^ in M" This article is devoted to another theory: the thin 
layer quantization. In this approach the particle moves between two 
equidistant infinite potential walls [7] or it is subject to some poten- 
tial force which in a proper limit makes it moving strictly along the 
surface [8]. 

The thin layer method seems to be natural for description of low 
dimensional motions in nanoelectronics. Recent progress in nanotech- 
nology caused a great activity in the field. Free particle energy spec- 
trum was investigated for thin layers around cylinders [9], tori [10] 
and arbitrary surfaces of revolution [11]. In this article we clarify the 
general properties of the thin layer method and establish its equiva- 
lence with Prokhorov quantization procedure [12,13] for codimension 
1 surfaces. For a vast majority of higher codimensional cases the thin 
layer method fails to yield meaningfull results [14], at least of the gen- 
eral form (1), while the Prokhorov quantization can be generalized 
directly as we show in section 6. 

In section 2 we describe the thin layer quantization method in a 
geometrically clear manner which allows us to deal with any codi- 

^The codimension is the difference between the bulk space dimension and the dimension 
of the submanifold. In particular, codimension 1 surfaces can be defined by one equation 
and have onc-dimcnsional normal space at each point. For these surfaces all notions of 
classical differential geometry are valid after obvious generalizations in the number of 
coordinates. For example, in the next section we would use the principal curvatures [6] of 
codimension 1 surfaces. At the same time, highcr-codimensional surfaces present certain 
difficulties, and we encounter some of them in the section 6. 
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mension 1 surface in M" without any complications in comparison 
with 2-dimensional surfaces in M^. In contrast with [8] we use only 
the most elementary notions of differential geometry. In section 3 we 
establish the equivalence of the thin layer approach with Prokhorov 
quantization method [12,13]. In section 4 a method of quantization 
with Hermitian momenta recently proposed by Encinosa [15] is anal- 
ysed. In sections 5 and 6 we consider surfaces of codimension greater 
than 1. 



2. Quantization for codimension 1 surfaces 

We consider (n — l)-dimensional smooth surface in M" and two 

infinite potential walls at the distance 6^0 from the surface. Free 
quantum particle moves in the thin layer of width 26 between these 
potential walls. We introduce a curvilinear coordinate system in which 
\xn\ equals the distance from the surface to the given point, and the 

coordinate lines of xi, . . . , .x„_i are orthogonal to coordinate lines of 
Xn- We have the boundary condition = ^'|^^^_^ = and 

Hamiltonian H = —^A with A being the Laplace operator in M", 

n n 
i=l k=l 

9ik-y^ ij 

where A^s is the Laplace-Beltrami operator on the surface Xn = 
const. One can prove [5] that g~^^^dng^^^ = divTz with "n being a 
unit normal vector to the surface, hence 

A = dn'^+div{lt)-dn + ALB. (2) 

Indeed, let us consider two surfaces, Xn = and Xn = e. Suppose we 

have an infinitesimal area dS at the surface x„ = 0. We denote the 
corresponding area element on = e surface by dS': 
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We have div{n) = + 0(e) = + 0{e), so g-^/^dng^^^ = 

div("n ). 

Now, let's take some point at the surface and consider another 

coordinate system yi, ■ ■ ■ ,yn- We choose it to be Cartesian and such 
that at the given point the tangent paraboloid of the surface is pre- 

n-l 

sented in its canonical form: Vn = ^ ^aVa- So, in the vicin- 

a=l 

ity of the chosen point (V = 0) the equation of surface is y„ = 

n-l 

^ ^02/0 + ^(2/0)5 ^a-s are the principal curvatures. The unit normal 

a=l 

is ria = , "f; + 0{yl) = kaVa + 0{yl), n„ = -1 + 0{yl) and 
'i+E kivl 

a=l 

n-l 

diYlt = J2ka + 0{ya). (3) 

a=l 



The surface be obtained by y — > y = y + en , 

and dy'^ = dya{l + eka + C(2/a))- It yields g 



n-l 

n (i+o(j/?))<i2/^ 

dS' _ a=l 

n-l 

n {l+Oiyi))dya 

a=l 



n-l 

n (1 + efca) + 0{ya) near the point 'y' = Q. At the line 2/a = Va = 

a=l 

1, . . . , 77, — 1 one has 

I n—l I /n—1 \ ^ n— 1 \ 

§ = i+^Efca+^^M (Em -E^^ (4) 

The relation (4) is valid at every point of the surface provided that 
one takes the principal curvatures at the same point. 
Following [7,8] we introduce a new wave function 



X{x) = -^{x] 

It is natural because 



dS' 

Is' 
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dV\^{x)\' = J dxnj dS\x{x)\\ 



\Xn\<S -S 
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so that the normal and tangential coordinates are completely sepa- 
rated. For the lowest energy solutions the nornal motion is restricted 
only to the factor of cos and the integration over Xn yields just 
the constant number. It means that the conservation of norm for x(a;) 
is satisfied. From (2)- (4) one gets 



+ 2dnX{x)dn + div "n • + div • x{x)dn 



For the lowest energy levels we have 

X{xi,...,Xn) = /(xi,...,Xn_i)cOS^ (5) 

and, after taking (5—^0 limit and subtracting an infinite (proportional 
to 1/(5^) energy, the Hamiltonian 

6 = -^A,s + ^i('j:K) -2Y^k^ (6) 



a=l 



is obtained. In simple cases the factorization (5) works very good even 
for not so small values of 6 [11]. 

What was presented before (infinite potential walls) is rather the 
approach of [7] than of [8]. But the difference is not very important. 
One could use an appropriate confining potential instead of infinite 
walls. It would lead to the lowest energy level function of the potential 
Kon/(x) instead of cos and to another infinite energy. 

Hamiltonian (6) contains quantum potential 




For two-dimensional surfaces in we get the result of da Costa [8]: 

R 



Vq = —^(ki — ^2)^- For a sphere ka = -k and the potential is Vq 



a^(n-l)(n-3) 
8R^ ■ 
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Physically this quantization can describe lower dimensional mo- 
tions in nanoelectronics provided that the restricting potential makes 
a layer of uniform effective width. Of course, it would be still a se- 
vere problem to check this potential experimentally due to the great 
energy^ proportional to the width of layer powered by —2. This energy 
is constant on the physical surface if the potential depends only on 
a distance from the surface. For more complex potentials the results 
may differ. One can get different infinite energies in different parts of 
the surface. It would mean infinite tangential forces, of course. But 
if variation of width becomes smaller and smaller while 5 — it's 
obviously possible to obtain some finite additional potential. 

3. Equivalence with Prokhorov quantization 

We should mention that there is one more method of quantiza- 
tion proposed by Prokhorov [12]. The motion of a particle is con- 
sidered as a system with two second class constraints but the only 
one condition is imposed on the physical sector: Pn^physix) = with 
Pn = —ihr^-^g^^^ . It means that 



Having solved some task by this method, one should put x„ = in 

the results after all the differentiations over ,t„ are performed. Due 
to (7) the probability to find a particle at the distance |a;„| from the 
surface does not depend on the value of x„, and we choose one value 
we need. For Prokhorov's view see [12,13]. 
From (7) and (4) we conclude that 



^Actually, it is not so great in comparison with the kinetic energy of chaotic motion at 

ordinary temperatures. It means that the method would not work well due to excitation 
of higher energy levels of the transverse motion. So, a cryogenic experiment is needed. 
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/ 1 / \ 1 

dl^physix) = %hys{x) I ^ ( 1^ 1 + 2 ^« + ^(^") 

Now (2) and (3) yield: 
- -^^"^Physix) 



So, the quantum potential coincides with the one obtained by the 
thin layer method. The reason is clear. The lowest energy level wave 
functions (in the model with two infinite potential walls) have nodes at 
Xn = ±5 and the bunch at .x„ = 0: 9„x = oi") equivalently, Pn'^ = 0. 
The methods of da Costa and Prokhorov are equivalent (disregarding 
the infinite energy of the thin layer quantization). 

Actually, one could use the calculation of A*(a;) from section 2 
and put dnX = (but note that in Prokhorov method one uses ^ as 
a wave function, not x)- 



= --^^LB'^phys{x) + 




1 \ n— 1 



Y^kA -2Y,kl\^phys{x). 



(1=1 / a=l 



4. On Hermitian momenta of Encinosa 

Recently Encinosa [15] proposed one more quantization method for 
a constrained free motion. The starting point is the Hamiltonian in 

fn-l 2 \ 

curvilinear coordinates Xj, H = \ \Y^ '^Pnj ' recipe is 

simple: pi — > pi = —ihg~^^'^dig^^^ followed by the thin layer method. 
According to these rules we have in our notations 



^ /(a;i,...,a;n-i)cos^ _ 1^ pj /(xi, . . . , Xn-i) cos 



7rx„ 
28 



^ hi{x) 

lh?f{xi,...,Xn-l) ra T^Xr^ 

at, cos 



25 



Effectively this result may be considered as zero quantum potential 
(geometric potential in the terminology of [15]). We should make two 
important comments on this point: 
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a) Quantization in curvilinear coordinates is dangerous, because 
the results usually depend on the choice of coordinate system. Nev- 
ertheless, the meaning of the method considered could be clear if 
momenta operators were self-adjoint, but it's not the case even for 
spherical coordinates on the sphere. 

b) Strictly speaking, the recipe is not defined correctly, because the 



operator odering problem in terms is not solved. It is not difficult 
to deduce the correct ordering for the zero potential theory from the 

relation -h^A = -/i^ £ £ g~^^''dif/^g"'dk = E r^'Hig^'^'g^'Pi in 

i=l k=\ i=l 

any orthogonal coordinate system. But this particular ordering is not 

natural a priori. Let's turn to the case of 5^ with stereographic coor- 
dinates xi = 2(R + X3) cot (I) cos (p and X2 = 2(R + X3) cot (|) sin ip. 
In [16] it was shown that in these coordinates the Laplace-Beltrami 
operator on the sphere of radius R + xs equals 



2 V 4(i? + X3)V 2 V 4(i? + X3) 

Even less natural it would seem for 5" in R""'"^: 



2 2V 4(i? + x„+i)2 



i=l 



xf . . . + .X^ \ ^ " \ / xj + . . . + x'^ ^"'^'^ 



4{R + Xn+i)y \ 4(i? + x„+i)2 



And, of course, for coordinates of section 5 in [16] the situation would 
not be better. 



5. Elementary cases of codimension > 1 

In general let's consider m-dimensional smooth surface in rep- 
resented by its tangent paraboloid at some point: 

Til TfX 

a=l 6=1 

a = m -|- 1, . . . ,n, = . In general n — m curvature forms 

k^"'^ can't be diagonalized simultaneously and the notion of principal 
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curvatures does not exist, but in this section we treat the simplest 
cases for which the diagonahzation can be performed. 

First of all, for 1-dimensional manifolds (curves) curvature forms 

k^"^ are just real numbers. In this case by a rotation in the space 
of ya one can get y2 = \kyl + 0(2/1), ^3, ■ ■ ■ ^Vn = 0{yf). The unit 
normal vectors are n^^^ = kyi + O^yf), = ~l + 0{yf), = . . . = 
ni^) = 0{yf); nf ^ = -6ia + 0{yl) for a > 3. Wc have li^'^^'^^'^ = 



6cti3 + C(yi) and after the transformation if if = if + enli 

a=2 



one gets dy[ = (1 + e2k + 0{yi))dyi, dy'^ = (1 + 0{yi))dya- 

Now we consider a smooth family of such coordinate systems and 
normals along the curve. We introduce a new curvilinear coordinate 

system in which xi is just the length along the curve while Xa = 
■ where is the minimal norm radius vector from the curve 
to a given point and li^"^ is taken at the same point on the curve as 
T^. In this coordinate system we have 

. _f{l + X2kf 
9ik-\^ I 

and A = Ac + A„ + [j^d2il + xsfc)) ^2 = A^ + A„ + j^d2 
where Ac is Laplace-Beltrami operator on a curve x^ = const and 

n 

An = ^ 9^ is Laplace operator in a hyperplane xi = const. 

a=2 

To proceed with the thin layer quantization we introduce a thin 



layer = 5"^ around the curve and a wave function x( 



X 



a=2 



VTTx^ *(x) such that / dxi J dSnormal\x{x)\'^ = I dV\'^{x)\'^: 

A^{x) = A^^M= = Aex(x) + Anx{x) + ^x{x) + 0{xa). 

After subtracting an infinite energy due to A^xix) it yields the quan- 
tum potential Vq = —^k'^ as in [8]. 

We should mention that it is very important for this result that 
the thin layer is spherical at every point of the curve. Indeed, suppose 
we have a straight line in M^. We can first embed a cylinder of any 
radius R into and after that, with much more thin layer, we embed 
a line into the cylinder. The first embedding results in = — 
and the second one changes nothing because the line is geodesic in 
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the cylinder. Note that it is possible to generalize our considerations 
to the case of embedding into non- Euclidean spaces in a sense that, 
if one has a free particle Hamiltonian (probably, with some quantum 
potential) in a Riemannian manifold, he may consider a thin layer of 
constant width around any codimension 1 submanifold. For a cylin- 
der it would literally reproduce the demonstrations of the section 2 
because in intrinsic geometry the cylinder is flat. 

As another simple example we consider 2-dimensional fiat torus 
isometrically embedded in 



p4. 



2 I 2 

2 I 2 
^3 + X4 



rI, 



We use the following coordinate system: 



f)i = arctan — , 

Xl 



(j)2 = arctan 



X4 



ri = ^Jxl + xl- Ri, 
r2 = \J xl + xl - R2 



with 



({n + Rif 






and A = At + A„ + 



— !— a 





{n+Ri? 









1/ 



:;:^^dr^ where A^ 



is the Laplace- 
const. We have 



Beltrami operator on the torus ri = const, r2 
the volume element dV = (ri -|- Ri){r2 + R2)dridr2d(j)\d(f)2-, the wave 
function x = y^(l + + ^ and the quantum potential Vq = 

-^{kf + k^) with ka 
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6. General theory 

In general the geometric construction is very similar to that of 
1-dimensional case. The unit normals to the surface (8) at "y^ = 
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are nf = <5«/3(-l + 0{yl)), rt^ = E k^^Vb + 0{yl), a,P = m + 

6=1 

1, . . . ,n, a,b = 1, . . . , m, and after the replacement 'y ~y' = 
y + Yl ^a'lT' one gets 

a=m+l 

m / n \ 

6=1 V a=m+l / 

dy'a = (1 + 0{yc))dya- It means that in curvihnear coordinate system 



gab{xi, . . . , Xn) = ^ I (5ca + ^ Xak^ 
c=l \ a=m+l 



ca 



n 

/3=m+l / 

with Xa = ~n^"^ • as in the previous section and fab{xi, . . . , x^) = 
gab{xi, ■ ■ ■ ,Xm,0, . . . ,0). In this system if some point has coordi- 
nates {xi, . . . ,Xn) then its position can be found by adding the vec- 

tor x^n to the radius- vector of the initial surface point 

a=ra+l 

with coordinates (xi, . . . If we choose Xa = ya at the surface 

Xa = 0, then at the Xa = hyperplane fab = Sab, x^ = (a and 

9db = Sab + '^Y^ ^ak^db + eaepkic^k^f^ . An easy computation yields 

a oc,l3,c 
oc,a a,P,a,b a,fi,a 

o,6 V a / a \ a / 



The problem is that the coordinates Xi are not orthogonal. Par- 
allel translations of the surface along one of the normals break its 

orthogonality to other normals. Indeed, we have = k^ + 

0{yc) and = 0{yc)- We suposc that all normals have unit 

length, hence = ^ • = E^i^yc^i? " ^ + ^(^c) and 

6,c 
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= Y^k^^k^^'uc + 0{y1). Now we see that parallel translations 
lead to violation of orthogonality condition because ■ = 

b,c ^ ' 

Wc introduce a new notation j^^} = jf^ = J2 [kfb ~ ^ib J ^bc^J/c 

6,c ^ ^ 

and find the metric tensor 



9ik 



I 



its determinant g = g + 0{e^) and the reciprocal tensor 



~ik 



f 



gab 

j:fLVg'% + 0{e' 

c,7 



c,7 



'c/3 



J 



Da Costa concluded [14] that the thin layer quantization would not 
work well in this situation because A contains terms with both deriva- 
tives da and da- 

The situation is quite different for the method proposed in [12]. 
All new terms in A have coefficients of order 0(e); the only suspicious 



term, J2i9ag''°')da 



I - E g^'f^^^Safs + Oie) ) da = 0{e)da, is 

not dangerous because /i"^ = 0. In Prokhorov quantization method 
the condition da [g^^'^^phys) =0) Va = m + 1, . . . , n is imposed. For 

a function x{^) ~ ij) ^{x) it means daX = Oj we get 

(f) (f) if) 



+ Xix) ( ^^«Vg) 9a ^ +0{Xa), 

a=m+l / U\ 



A„ = X] dl- It is easy to see that -^da^/g = -2g^''^dag'~^''^, and 

a=m+l 
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the Hamiltonian is H 




a=m+l 



n 




=0 



One more (not too difficult) calculation yields 




(9) 



The thin layer method would not give this answer because in this 
method we have d^X — only at the original surface and generally 
daX ^ ^: so that terms with daX are not negligible and factorization 
analogues to (5) is not a good approximation of the exact thin layer 
solution. It's not surprising. Let us take a small element of the surface 
and its (^-neighbourhood, 5^0. The well-known theorem states that 
J dVA'^ = J dS ■ x/"^. Wc have A^I* ^ f dS^V^ ^^^^ -^^ codimension 
1 case the normal projection of leads to p term in (because 

~ I and V ^ S). Tangential components of V* result in fi- 
nite values of A* due to dS± ~ S. In general case V ~ §n-m ^^^^ 
dS ~ (5"~"*, hence the tangential components of yield finite terms 
in A^' again; but normal projections gain some components orthogo- 
nal to transverse hyperplanes (xa = const). The corresponding angles 
are of order S but ~ |. It results in finite terms in A^'. So, 

normal components of V* infiuence the tangential dynamics. Hence, 

the thin layer method of [7,8,14] does not yield the result of the gen- 
eral form (1). It means that the Prokhorov method in such cases does 
not correspond to the motion in the uniform thin layer but it turns 
to be more powerfull as an abstract quantization method for second 
class constrained systems. One could find, of course, some orthogonal 
coordinate system in the whole vicinity of the surface with normals 
depending on Xa and, may be, he would be able to determine a ge- 
ometry of the thin layer for which some factorization analogues to (5) 
would be a good approximation. Quantum potential in this approach 
is likely to coincide with (9) at least with a certain realization of it, 
but generally the layer would not have constant width and the set-up 
of the quantization would not be as easy and clear as the original one. 
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